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ABSTRACT 



In a recent paper* by T. J. Harvey, there is presented the idea of 
employing known solutions of free vibrations of nonlinear single-degree 

conservative systems to obtain solutions of forced vibrations of these 
systems. The forcing function is taken to be of the same form as the 
free response and the solutions are obtained quite easily. Harvey makes 
the point that while this excitation is not the sinusoidal excitation 
that many analysts have attempted to study, for many purposes there is 
reason to prefer some alternate form of excitation if by doing so useful 
solutions are readily obtained. In a paper of several years back, M. A. 
Sadowsky* derived the characteristics of several nonlinear spring systems 
for which the free vibrations could be obtained in rather simple closed 
forms. In particular, Sadowsky introduced two new springs which give solu- 
tions in terms of the elementary functions and several others which give 
solutions in terms of the Jacobi elliptic functions. 

Since there are actually only a few fully developed cases of forced 
nonlinear vibrations available to analysts for testing the power and ac- 
curacy of analytical methods and for other purposes, it appeared to be of 
interest to apply the viewpoint of Harvey to several of the particular 
cases given by Sadowsky so as to work out all the details of the corres- 
ponding cases of forced vibration and thus afford several new examples. 

The major content of this thesis is the full development of these cases. 

In addition, however, the viewpoint of Harvey is applied to two cases 
of piecewise linear springs simply for the purpose of illustrating such 



*See references on page 41. 
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applications and for the purpose, as well, of obtaining additional cases 
of nonlinear vibrations. 

Also, since the forcing functions differ from the sinusoidal func- 
tions that many analysts have used, it was considered desirable to obtain 
and exhibit a simple measure of the degree of "distortion” of the actual 
forcing function as compared to the sinusoidal shape. Likewise, the dis- 
placements are not sinusoidal functions of time, and the "distortion" of 
the response (as compared to a sine curve) is also examined. 

This work was done largely during the spring of 1961 at the U. S. 
Naval Postgraduate School, Monterey, California. The writer wishes to 
express his appreciation to Professor J. E. Brock for having suggested 
the topic and for guidance during its development. 
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SYMBOLS 



t * time s independent variable, chosen so that y(0) = 0 

Q X- s dimensionless independent variables 

^(r) s deflection, normalized 

p(f) * external force applied to the system 

Si tf) = spring function such thatS(o) = 0 and Sfa) - 1.0 
$(#) has the same dimension as does q- 

= assignable dimensionless constant as defined in the relation (5) 

■)fl = mass or moment of inertia of the system 

s spring constant for a small deflection 

l/Jf, = \R/n\ = natural frequency for a small free oscillation 

Tv s ~ P er i°^ a small free oscillation 

To s period of a free vibration ( p = 0) 

U = frequency of vibrations 

J * s period of forced vibrations 

ft -$(%) - amplitude of deflection 

— denotes the evaluation at "£ - % 

c( - modulus of Jacobi elliptic func* ion when a solution contains 
such functions (refer to Table 1) 

oC.j - modulus of Jacobi elliptic functions whe^ a spring function 
contains such functions (refer to Table 1) 

K = the complete elliptic integral of the first kind with modulus <Y 

Kj - the complete elliptic integral of the first kind with modulus % 

/\ b the distortion of $t) wave form as defined in Eq. (17) 

B = the distortion of F£t) wave form as defined in Eq. (18) 

(9 = a constant for the springs of cubic parabola types such as 

RlYj - Y+P Y J 

Y s deflection which is not normalized 



V 



LIST OF ILLUSTRATIONS 



Table 

1. Solutions for ten different nonlinear springs adapted 
from Sadowsky (Ref, 1). 

Figure 

1. Approximate spring characteristics adapted from 
Sadowsky (Re f . 1 ) . 

2. and 3, Explanations about measures A and B, respectively 

4. Frequency response curves, spring characteristic, and measures 15 
A and B for S(})= 

5. (a) and (b). Wave forms of %$) and f &) » respectively for 16 

6. Frequency respon°e curves, spring characteristic, and 
measures A and B for S(‘J) = 

7. (a) and (b) . Wave forms of tf(t) and fit) , respectively for 

S&) = MyuJfy 

8. Frequency response curves, spring characteristics, and 24 

measure A and B for S(fy) = & + 'k if? 



9. (a) and (b) e Wave form of iffc) and p(t) , respectively for 25 

sm = %+kf 

10. Frequency response curves, spring characteristic, and 
measures A and B for 

11. (a) and (b). Wave forms of and f(t) , respectively 

for Sft) = 

12. Frequency response curves, spring characteristic, and 34 

measures A and B for a piecewise linear softening spring. 

13. (a) and (b). Wave forms of %(t) and f(f) , respectively for 35 

the spring in Fig. 12. 

14. Frequency response curves, spring characteristic, and 39 

measures A and B for a piecewise linear hardening spring. 

15. (a) and (b). Wave forms of $(t) and f(t) , respectively for 

the spring in Fig. 14. 



1. Introduction,, 



The equation W s + A $ ~ F c o"< u’T (I) 

describes the behavior of a linear single-degree-of- freedom conservative 
system subjected to a sinusoidal excitation, 

A natural extension for the study of a nonlinear conservative system 
is 

n £ + k s &) = h tit ( 2 > 

where $($) Is a nonlinear spring function, such as -S^) - ^ % . Many 

investigators have studied Eq. (2) in detail and obtained solutions which 
indicate the effects of parameter CO upon ^ . But Eq, (2) is complicated 

to solve analytically, and in most cases the solutions are obtained by 
approximation methods. 

T. J„ Harvey (ref. 1) suggested a simple analytical method of study- 
ing the forced vibrations of nonlinear systems by making use of the solu- 
tions for free vibrations, i.e,, where the right hand side of Eq. (2) is 
zero. Some years ago M. A. Sadowsky (ref. 2) presented simple solutions 
of free vibrations for various nonlinear springs. 

This thesis applies Harvey's method to several of the cases given by 
Sadowsky, thus obtaining several interesting solutions for forced vibra- 
tion of nonlinear systems. 

In addition, a study is done on the following: 

(a) The wave form of %(£) is compared with sinusoidal one which 
is that of a linear system, 

(b) Traditionally the right hand side of Eq. (2) is a sinusoidal 
function, as shown. By Harvey's method it differs from a sinusoidal func- 
tion, but is still a periodic function. These differences are compared. 
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(c) Harvey's method can be applied to other than Sadowsky's 
nonlinear springs. It is applied to the spring characteristics which 
consist of linear segments, as examples. 
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2 c b^sic equations. 

The forced vibration of an undamped nonlinear single-degree-of- free- 
dom system is described by the equation 

'fot % + i$Cf) ■■= p(0 (3) 

where a nonlinear spring function normalized so that $(C) = 0 

and $'{&)** hO , and f(t) 1 s the external force applied to the system* 

Eq. (3) is simplified by substituting: 

and i 





< 


Then 


+ 


Eq. (3a) 


as it stands 



(4) 

(3a) 



periodic and if we presume (by way of analog^with linear systems) that 
slight damping present in the system but not accounted for in the mathe- 
matical description causes the transient" to die out, the "steady state" 
solution is generally not obtainable except by difficult analyses leading 
to approximations. However, the analysis is greatly simplified by henceforth 



(3) 



confining attention to only such f(i) as satisfy the relation 

* 

I % 

where is a dimensionless constant, either positive or negative. Such 
a forcing function is called a "natural forcing function" (ref. 1). A 
positive means the external force is in-phase with restoring force, a 

o 

negative means out-of-phase (180), and f>- 0 means free vibration. 

Substituting the relation (5) into Eq. (3a), we obtain 



■jfr = o 

This is further simplified by putting 

a — x. 



( 6 ) 



(?) 
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Then 



( 8 ) 



M 

ile* 



+- Jis 



which is described as the case of the pseudo- free oscillation. To solve 
Eq. (8), let the initial conditions be 

*•» = ° • M = K 

The potential energy stored in the spring is 

u(?)= [*$(?)& , u(o) 

SO 

Furthermore, by conservation of energy, Vc is related as 

v; a = 2 uQ) 

where ft is the amplitude of deflection. If T/4 is defined as the quarter- 
period of oscillation, ft can be written as ft - ft( } 

Eq. (8) is integrated and simplified as 



(tI) ^^“(b-iacft) 



(9) 



( 10 ) 



Integration of Eq. (9) leads to 

ft dft _ 

e - l 

Again Eq. (10) can be written in the following general form with r 
as a parameter, 

B — &C ft J ( 11 ) 

The inverse form of Eq* (11) can be expressed as 

Since the dimensionless variable Q is closely related to time t, 
using Eqs. (4) and (7), 

0 - I- f> ( 12 ) 

Then $ = 1 ) (13 > 

This is the solution of the forced vibration of Eq. (3) in which the 

external forces f(t) satisfy the relation (5). For free vibrations 



and 






Then the solutions of free vibrations are from Eq. (13)< 



? “ ? ( «vr. f ) 



(14) 
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Comparing Eqs, (13) and (14), the solutions of forced vibrations are 
obtained from those of free vibrations by multiplying the independent vari- 
able t(^ time) by the factor \j 7 — , The assignable dimensionless con- 

stant p is the ratio of the external force to the restoring force, namely 

F 

kStf) 

where the relation (5) is rewritten. The amplitude of external force p 
is given as 2.0^ or 3,5^ , for example, which means the force unit corres- 
ponding to the value of the spring constant (at small deflection) multiplied 
by 2,0 or 3,5 units of deflection respectively. Or in other terms, Eq, (3a) 
represents a forced vibration wherein plays the role of a restoring 
"force" and ^ plays the role of an exciting force; the time scale 

has been chosen so that both "forces" have the same dimension as the "de- 
flection response", J f(t). Thus we can say that the amplitude of exciting 
force ( F/lk ) has magnitude 2,0 or 3,5, say. For every amplitude of de- 
flection f , the values of o — 'S(tf) are found, which, in turn, deter- 
mines the assignable constant ^ . 

Sadowsky's paper gives the essential details for ten cases of nonlinear 
springs for which the free oscillations are relatively simple. Two of the 
solutions are in terms of elementary functions, and the others involve the 
Jacobi elliptic functions. Corresponding to each of these cases, by using 
Harvey's method, it is possible to construct a complete solution for the 
forced oscillations resulting from "natural" excitation. In this thesis 
four of the cases are developed in detail- -the two cases involving only 
elementary functions and two of the cases (cubically hardening and cubi- 
cally softening -- i,e., Duffing* s spring) involving elliptic functions. 
Also, two cases of piecewise-linear springs (cases not included among those 
given by Sadowsky) are also developed. 
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Sadowsky v s solutions are expressed in terms of the initial velocity 
\/ c , which can be replaced by V , for convenience, with the energy rela- 
tion V/ — * in the conservative systems. The solutions are tabu- 

lated in Table 1 and their approximate spring characteristics are shown 
in Fig* lo The symbol of represents the modulus of Jacobi functions, K 
means the complete elliptic integral of the first kind with modulus c< 

In Fig. 1, the numbers in circles represent the spring numbers as in- 
dicated in Table\ Springs marked as 1, 2, 3, and 4 are drawn heavily. 

The family springs 7 is in between springs 1 and 4 depending upon the modu- 
lus o( i , where for — <2 it coincides with spring 1 and for 
it coincides with spring 4. The space between springs 4 and 3 is filled 
with the family of springs 6 depending upon the modulus Of; . For o {\~Q 
it coincides with spring 3, and for C(r-l'0 it coincides with spring 4. 

Each one reaches a maximum and drops down to zero. The family of springs 
5 is of softening type, and so is the elementary spring 2 obtainable with 
^1*0 . Springs 8, 9, and 10 are cubic parabolas, they are not shown 

in Fig* 1 since they can readily be visualized. 
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Table 1 - The solutions of forced vibrations for Sadowskys springs 
(adapted from Sadowsky; ref. 2) 

The solutions of forced vibrations such as 

VI $ + k 5 (t) - Ffc) 

are tabulated below, in which f~(t ) satisfies the relation (5). 
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Fig. 1 Characteristics of Sadowsky's springs, 
(from Sadowsky, ref. 2). 





Fig. 2 Distortion of the 
wave form. A is 
defined in Eq. (17). 



Fig. 3 Distortion of the 
wave form. B is 
defined in Eq. (18). 
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A, Frequency response. 

When the exciting force is such that the relation (5) is satisfied, 
the period of forced vibration is readily obtained from the solution of 
free vibration. If the period of free vibration is obtained as: 

T c = 



then by comparing Eqs. (13) and (14), we see that 



T 

To 



I 



'l I - p 

Since l 0 is varying with ^ , the ratio 

period of a small free oscillation Tit 



(15) 

may be related to the 



_T 

Tn 



* \tJ 



i 



n-> <l6 > 

This means that corresponding to a single frequency response curve 
of if versus , obtained for free vibrations, we obtain a full 

family of response curves for forced vibrations, i.e,, curves of ^versus 
^-1* ^ for various values of ^ 
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B. Wave form of l ^(t) . 

When the exciting force is such that the relation (5) is satisfied, 
the solution is exactly the same as that for free vibration except for 
the change in frequency, which is times that of the free vibra- 

r 

tion. Therefore the wave form of depends upon the amplitude of 

deflection only, and not upon the amplitude of exciting force, When 
is small, the springs are very nearly linear, and the wave forms 
are almost purely sinusoidal. As $ increases, the wave forms of $(t) 
are distorted from the sinusoidal forms. 

Furthermore when the spring characteristics are symmetrical, the 
wave forms are quarter-wave symmetrical. In addition, they tend to lie be- 
tween a rectangular and a triangular wave form. 

To show the distortion of ( | )(t ) wave forms, a measure A is defined as 

X - ■ f * 

If we let the quarter period T/4 correspond to % , the measure A 

represents an area under the curve of M as shown in Fig. 2. Note 

% 

that A— / » for a linear spring or a small oscillation. 
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Co Wave forms of p (ft) . 



When the external force f-(t) satisfied the relation (5), the wave 



cal in the same manner as that of tyft) . For a small oscillation fft) 
is almost a pure sinusoidal form, from which the wave form is distorted 
as ^ increases. 

To compare this wave form with the traditional forcing function of 
sinusoidal form, a measure B is defined as 



form of \~ft) is the same as that of the spring force Therefore 

it depends on ; only, not the amplitude of the external force p 

For a symmetrical spring, the wave form of F ft) is quarter-wave symmetri- 




Therefore 




(18) 



Note that: 3— /'£> for a linear spring or a small oscillation. 
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3. Examples. 

The results of six examples are presented in Fig, (4) through Fig, (15) < 
Figures 4, 6, 8, 10, 12, and 14 are all similar representations as follows: 

(a) Ordinates are ' (or where C is a constant). 

(b) Solid curves represent frequency response curves for 

labled values of p- > including p - o > the case 

of free vibration. 

In-phase forced vibrations are indicated by positive values of p , 
out-of-phase by negative values. These curves relate (ordinate) to 

(bottom scale of abscissa). 

(c) Distortion of is indicated by a measure A, distortion 

of p(t) by a measure B, both functions of ' alone, and 
plotted on vertical scales. 

(d) The dotted curves show spring characteristics S(if) (top 

scale of abscissas) as functions of ^ . In Figs. 12 

and 14, this relation is shown to smaller scale in the 
inset and not in full scale as in the preceding figures. 

Figures 5, 7, 9, 11, 13, and 15 actually display the distorted wave 
forms for ff(t) and f(t) . 
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Ac Example 1 






.■ini j- 
U-J. * u 



The characteristic of the spring is extremely hardening as shown in 
Fig. 4 (dotted curve). The solution of the forced vibration is already 
tabulated in Table 1, namely 

Ami ‘f =- Orf- uit ( 19 ) 



where 



U) = « A* f) p-~/° , f = -£- 

(a) Frequency response. 



( 20 ) 



From Eq. (20) we can get the solution which indicates the 






effects of frequency upon if- , making use of T~ io and Tn — 



_ -ZT 
Ufa 



In 



the dimensionless form, we get 



-X- — coi 

Tn 






e 



( 21 ) 



For the given value of p , f 3 is obtained by dividing 
p by S(^J > then the ratio is obtained from Eq. (21). The fre- 

quency response curves are plotted in Fig. 4. 

(b) Wave form of tf(T) . 



( 22 ) 



Eq. (19) is rewritten in a convenient form, 

%(t) _ 4am~' A^n iv t) 

For a small oscillation, = ***! U't 

r*° $ 

With Eq. (22) the wave forms of fr(t) are plotted for ^ * 0, 
0.5, 1.0, 1.30, and in Fig. 5(a). Since Eq, (22) shows that %(t) is 
quarter-wave symmetrical, plotting is done for a quarter period. For a 
small oscillation, it is clearly a pure sinusoidal form, and as in- 

creases it is distorted. To show the distortion, the measure, A(as defined 
in Eq. 17) is found. The integral is not easily performed analytically 
so that a numerical integration method was used. The measure A is a 
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function of ^ , and is plotted along ^ in Fig. 4. 

(c) Wave form of f?t) 

Using Eq. (19) and the relation (5), 

h (*) «= f>kS [$(V] tod f - f> k S ( f ) 

rearranging these into a dimensionless form, we get 

pit) ^ i urt 

f 0+ ^ 2 } u ^ viy* ( 23 ) 

(j > Fit) > 

For a small oscillation, fzF — ( ° 6 

ft F 

which is purely sinusoidal. As i increases, it is distorted. With Eq. 
(23) the wave form of p(t) is plotted in Fig. 5(b). 

To compare f(jt) wave form with the traditional sinusoidal function, 
the measure B(as defined in Eq. 18) is obtained: 

which is a function of ft only, and is plotted in Fig. 4. 

(d) Discussion. 

This spring has an extremely hardening characteristic. At 
ft — ^ , the spring force becomes infinity, its slope is infinity, and 

the measure B is zero. For slightly smaller values of ft , the met sure 
B is far away from unity. Fig. 5(b) shows that the wave form of f(t) may 
become quite distorted. 

When ft is far away from the value of % , the measure B is not far 
away from unity, and Fig. 5(b) shows the wave form of f(t) is not so much 
distorted. 
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B. Example 2< 



Sfif) = t 
0(3 > u>ih 3 it 



This spring has a softening characteristic as shown in Fig. 6 0 

The solution of the forced vibration is already tabulated in Table 1. 

rj j== “- sXtsfiii AsCfi k T ^ 24 ) 



where 



U) (u.l-i 4*i f JJ l-p 

(a) Frequency response. 



k S 



From Eq. (25) 



x = «4?=k_ 

I* 'i I— p 



(25) 



(26) 



For a given value of p , p is obtained as a function of 
~jj_ , and the dimensionless values of ~pj~ n are found. The frequency re- 
sponse curves are plotted in Fig. 6 £orf=a and db&lfe. 

(b) Wave for of fl(t) . 



(27) 



Eq. (24) is rewritten in a convenient form, 

cffrj -dirtJl ft Akt ‘Jt’) 

& 0 . Utf) , , 

For small oscillation, 'ift’X-c <kT 

which is purely sinusoidal. With Eq. (27) the wave form of ft(t) is 
plotted in Fig. 7(a) for , 1.0, 2.0, 5.0, and 00. As ft increases 

it becomes distorted. And as ft approaches infinity it goes to a rec- 
tangular form. 

To show the distortion, the measure A is found and plotted along ft 
in Fig. 6. The integration to find the measure A is performed by a 
numerical method. 

(c) Wave form of J-(ir) . 

Using Eq. (24) and the relation (5), we get 

Fin = pks[it(»l , F-pksn) 



rearranging these into a dimensionless form, 
f(t) u/t 



p ( I - u.ty /j - 



(28) 
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For a small oscillation, 




F 



which is purely sinusoidal. As if increases it is distorted. With Eq. 
(28) the wave form is plotted for = 0, 0.5, 0.66, and 1.0 in Fig. 7(b). 

To compare f(t) wave form with a sinusoidal one, the measure B 

is found. 



its maximum value (which is equal to 0,384). Beyond this point it decreases 
and the spring force approaches zero as < approaches infinity as shown in 
Fig. 6. Up to - 0.66, the measure B is not far away from unity, and 
Fig. 7(b) shows that the wave form of F(t) is not distorted appreciably. 

Far beyond the point ^ = 0.66, the value of B is much greater than 
unity and the wave form is greatly distorted. 




(29) 



which is plotted along if in Fig. 6 
(d) Discussion. 



The spring is of softening type: at if * 0.66 5$) reaches 



18 



C. Example 3. Duffing's equation. Sift) “ (f ' r ^ d 

When the spring characteristic is a cubic parabola, the equation of 
forced vibration is 



*r + -ft(Y HX S ) = Fit) ( 30 ) 

in which Eq« (3) is rewritten here. X the deflection. If it is ex- 
pressed in inches* p> is in inch and in lbs/in. 

To solve Eq. (30) with the right-hand side set equal to zero* the 
general cubic parabola (Y+ pX j can he expressed in the form of Table 1* 
l.e.* ( + by proper substitution regardless of (3 0 

Lst T « 31 > 

Then the equation of free vibration becomes 



jj- + -~o 

in which ft is a dimensionless quantity. 

sir) - i t ( 32 ) 

Then we obtain the solution of forced vibration immediately from 
Table 1. 



where 



o< 



'ft(t) - } On (u.'t- K , J 

uj 



y 

°< = }(< ) as the modulus of Jacobi functions 



(33) 



(34) 



p JJjq in a dimensionless form 

And fC is the complete elliptic integral of the first kind with 



i.e. 



(% 



K = ( i- AwQ) * 

>c y 



(a) Frequency response 

4-K 

From Eqs. (33) and (34) with -g- 

T _ , 1 /£ K \ -=L 

T ~ J I + ^ ' it J J i~ 



we get 



r 



( 35 ) 



In Fig. 8 the frequency response curves are plotted for p - o s i 



and ± 2 . 
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(b) Wave form of ^(Z) « 

Eq. (33) is rewritten; 

m 






- 6, ( ( ict- K , 0(J 



For a small oscillation. 



Lt. //i 






t 



ML 

f c i 

which is purely sinusoidal. As if- increases it becomes distorted. 



However there is a limit in the distortion. As $ approaches infinity the 

i 

value of cK does not exceed fj- . In Fig. 9(a) the wave form of if(t) 
is plotted for § ** 0 and ft = Cxo . 

To show the distortion the measure A is found 9 A 
which is plotted in Fig. 8 along ^ . It does not change appreciably 
from unity. Even if ft approaches infinity, A approaches 0.942. In other 
words, the wave form of ft(t) is not distorted much from the pure sinusoidal 
form. 

(c) Wave form of p(t) • 

From Eq. (33) and the relation (5), 

F(t> « fj* S f m) ' f = PrTf s ( V 

rearranging these into a dimensionless form, we get 

-&P- - ... 1 _ ! cnfait-K) + if 

F < 

with the modulus o< 

n. ML . ^ 

Jc±m = - jMV± 



(36) 



For a small oscillation. 



which is purely sinu- 



soidal. With Eq. (36) the wave form of f(t) is plotted for ft «* 0, 1.0, 
2.0, and infinity in Fig. 9(b). 

To compare p(t) with a sinusoidal force, the measure B is found. 



p _ _ZL ( £±± V* 
v ~ a K I 4 ) 



which is plotted along ft in Fig. 8. 

(d) Discussion. 

This spring is a hardening type as shown in Fig. 8 (dotted 
curve). Every type of cubic parabola, such as Y can be reduced to 
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the form Slfr) — + 2 <f by letting Y ^ without losing the 

generality of the governing equation (Eq. 3). 

Fig. 9(b) shows that the wave form of f (t) is not distorted 
to an appreciable amount, and the measure B is not far away from unity. 

The application of Harvey's method to Duffing* s equation 
was considered in a recent note (ref. 3) which appeared after this thesis 
was nearly completed. For the hardening spring, considered in this section, 
the results given in the note agree with those given here. However, in the 
case of the softening spring, treated in the next section following, the re- 
sults do not agree since the note contains an error. In place of the cn 
function in the note, the cd function should appear for the softening case, 
and the argument of the function should also be slightly modified. 
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Dc Example 4, Duffing 6 s equation, ft - 5 - ft* 

When the spring characteristic is a softening type of cubic parabola, 
the general characteristic equation is expressed as 

Rm = r- fir 3 

which can be reduced the form of Table 1 by letting Y = ft , where 

ft is a dimensionless quantity. Then without losing the generality of 
the governing equation (Eq, 3), we get 

sm = 



The solution of forced vibration is already tabulated in Table 1, 

ft(t) = J M(uXM) 



(37) 



where 



(a) Frequency response. 



(38) 



~-±K. 



From Eq. (37) and (38) with f = ~u^~> we § et 



Tn 






(39) 



where K is the complete elliptic integral of the first kind 0 The fre- 
quency response curves are shown in Fig* 10 for p s 0 and ± 0.2 
(b) Wave form of ft(f) * 

Eq, (37) is rewritten in a convenient form 2 
__ M(lctM) where o( = ft & 

For a small oscillation, 

which is purely sinusoidal. As ft increases it is distorted as 
shown in Fig, 11(a), 



To show the distortion the measure A is founds 

A -(*)40+iSrJ* 



(40) 
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and it plotted in Fig. 10 along *;> 0 

(c) Wave form of f(t) . 

From Eq. (37) and the relation (5), 

R0 - 5 M and F-fjif SC r t 

rearranging these into a dimensionless form, we get 



Eft) 



- | _ °0 — (titM'A 



(41) 



For a small oscillation, ^ ~ 

which is sinusoidal. As increases it is distorted as shown in Fig. 



11(b), 



To compare this wave form with a sinusoidal one, the mea- 



sure B is found, 

6-(£)Ad+ + 




(42) 



which is plotted along £ in Fig. 10. 

(d) Discussion. 

This spring is of softening type, and a t ^ * 0.82 it reaches 

its maximum. After this point it drops down, and at the spring 

force becomes zero. 

The measure B does not change appreciably up to 0.82, 
and Fig. 11(b) shows that the wave form is not distorted much. 

For any value of p which indicates the nonlinearity of the 
spring, the above solution can be applied. 
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E o Example 5o Piecewise- linear softening spring** 
Suppose a symmetrical spring 
(as shown) is described by the equations, 

(f) — kt (f |</l^ (- 

ip c 



S(» 




> 



5(f) - h f + ( k -fu)c 

and it is considered the case of k^Oouly 
The equation of free vibration is 

m# + sit) c^3> 

When the amplitude of deflection ft is less than C s it is linear, 
and solution will not be included here* The case ofjf^C will be consider- 



ed* 



Let us denote for 


o 4 1 4 t c 




f ° r 


ilc - 1 4 ^ 




where tc. - the time when Vfe) ■ C 

and let fyn\ and UjF = 






Then Eq 0 (43) becomes 






i + “'■» -o 


o 4 t i 




£ + n£&*(*?-ui)c -o 


t C ^U T 4 


(43a) 


The general solutions of Eqs, (43a) are 


* 




ut)= n 


0 4t4 *c 






IF*} 


(44) 



* In section E and F (which follows), for convenience in notation, 

only the first quarter cycle of motion is considered* 
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where Q, J) z , <f> , and (p z , are arbitrary constants, which together with 
tc may be determined by the following boundary conditions; 

ty(c) — C as an initial condition. 

ffio) ~ -2 ^(i) by the conservation of energy. 



Solving for the constants in Eqs. (44), we get 

n = + Sa-cj f 

a=-cci--g; + l («) 

4 — -b, 

<f>, - 

To obtain the frequency response curves, the period of free oscil- 
lation ( ^ “ 0) To can founds 

kx)“> 

because ^ 2 (t) reaches its maximum value of ^ at t- „ Solving for J 0 

S«: T. = |r (J - <p 2 ) 



we 



where T n = ~r is the period of a small free oscillation 

By making use of Eq. (16) for forced vibrations, the period of forced 
vibration is obtained, 

/ id \ / , $1 A I 

(46) 






* - (*K 






■n^ / J »-/> 

As a specific example let us take ki — 4- kz in the above symmetrical 



spring. Then we get U), = 2u £ , and the solutions of Eqs. (44) become 

0 4t^tc 

&(t) =(f+3c) Al*(aut + h) —.3 c 
X{= [c x + 2c(t-c) + ~(2-cfJ y * 



(44a) 



where 
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'U <J, ( P; ) 

^ = /M'r-ffe) - k 

(a) Frequency response. 

From Eq. (46), 

_TL 

r* ( r ; 



(46a> 



>/ I — f> 

where is an assignable constant as defined in the relation (5). In 

Fig. 12 the frequency response curves are shown for p=r o , ± o.$ Jz,C » 

and i k,(L . 

(b) Wave form of ^ f ) . 

Eq. (44a) is rewritten. 

= J- 4^1 U\T o 4 1 ^ t t . 

mi -[>+ - f tn <“•>> 

For the motion is purely sinusoidal* and as ^ in- 

creases it is distorted* At around ^=-Z*jrc , distortion reaches a maxi- 
mum* and for further increase in £ the wave form approaches a sinu- 
soidal form again* In Fig. 13(a) these two wave forms are shown: one is 

sinusoidal and the other shows the maximum distortion at „ 

To show the distortion the measure A is found. 



A = j ( i - yrrj r + (<+¥) 






(47) 



which is plotted in Fig. 12 along % /c • Actually the distortion is so 
small that the measure A is not much greater than unity s where the maxi- 
mum measure A = 1.023 at ft ~ 2 

(c) Wave form of F (t) 

From Eq. (45) and the relation (5), and with the spring 
force 5(%-) = {3C+tf) , | £b c , for this spring, we get 
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( 48 ) 



F(t> 4 ft 



F 
fit) 



din U), f c 



c ^ f < r 0 



= A*i (ti/t + § 2 ) 

In Fig, 13(b) the wave form is plotted for > ^ 

and i'c • Since the spring force has a slope discontinuity at [f-c » 
the wave form of f-(t) is slope discontinuous at T^f c , 

For it is a pure sinusoidal form. As !f- increases 

it becomes distorted; maximum distortion occurs at —2$c approximately. 
Further increases in “ bring out a sinusoidal form again, similar to 
}f( t) wave forms. Every wave form of f(t) does not go beyond a limiting 
curve which is marked as "Envelope" in Fig. 13(b), 

To compare the wave form of f(i) with a sinusoidal form, 
the measure B is found, 



which is plotted along ^ in Fig. 12. At ]f-— 2 *5C , approximately, j5 

reaches its maximum value of 1.22. At this point the distortion is a 
maximum. 

(d) Discussion. 

The solution of free vibration is found by a classical 
method, and for the forced vibration Harvey's method is applied to that 
of free vibration. 

Fig. 13(b) shows that f(t) is not far away from a pure 
sinusoidal form, and the measure B is not much greater than unity. 
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F. Example 6« Piecewise linear hardening springe 



Let is s consider a symmetrical 



For free vibrations. 




where 



»> t + sit) - ° 

Let the initial conditions be 

$(£) , fro) = * u( f) 

Then the solution is obtained from Eqs, (44) and (45) ? 

%(t) - d, <Lui u),t a i* tit Tc 

tit) — l • 4 '-) &*( rfyt.) t \c tet i 4 % 

D, = [c i + ^ccf-cj+4r|^;j 



(50) 






~^ c ~~ 4', j 

^ =W fe|7x) -***"&) 

(a) Frequency response. 



we get 



(51) 



F - 0, 

c Vt. 



From Eq. (46) and with 2ii\=uJ^ 

T_ _ (± _ iL ) — 1— . 

T» ( ^ ' »l i - /© ' 

In Fig. 14 the frequency response curves are plotted for 
i 4<C 9 and ±S$iC' For » it is linear and the ratio 

) is unity. And as ^ becomes much greater than C „ the spring 
characteristic again approaches linearity and the ratio ( ) approaches 

0.5. 

(b) Wave form of %(t) . 

Eqs, (50) are rewritten: 



1£> = JL ^ ^ t 
t t 



(52) 



-*y- - (I” X 



In Figo 15(a) the wave form of %(£) is shown for ^^4 
and jf—iyc* For % the spring is linear and the wave form is purely 
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sinusoidaic When ^ becomes greater than C , the wave form is distorted, 
at approximately » bffc the distortion reaches a maximum,, For farther 

increase in y , the spring characteristic approaches linearity again* 
the wave form approaches a sinusoidal form 

To show the distortion the measure A is found. 



A-wn(- $)('-“*«*■) + •#$(>>%■) 

(%)(<- #<r) < 53 ' 



1, 



which is plotted along f/f in Fig. 14. At ft ~i s'C the measure A reaches 
a minimum (A = 0.965). 

(c) Wave form of f(t) . 

The spring characteristic is 

0['f) - kzf-^kiC for 



tiii d.t 



0±t it 



Substituting Eqs. (50) into the above equation we 
P/t) = % 

F - ic 

F0 - ^(44 /><£.) 14 

F~ 



( 54 ) 



In Fig. 15(b) the wave form of f(r) is shown for 8 

i-.S'C 8 and 3.0 C . There is a limits beyond which the wave form 
does rot go, as indicated as "Envelope" in frhe figure. 

To compare the fit) wave form with a sinusoidal form the 
measure B is found. 



(55, 

which is plotted in Fig 0 14 along « For the spring is linear and 

J~(t) is a pure sinusoidal form, so the measure B is unity 0 As goes 

beyond C 9 the measure B decreases rapidly, and at ^ *l»5c 9 approximately, 

B s 0,735 as a minimum* For further increase in ^ 9 the measure B 

slowly approaches unity where the spring characteristic becomes linear ageim 
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(d) Discussion* 

This example is similar to the previous one* from 'which 
the general solution is obtained* The measure B and the wave form of 
F(t) shown in Fig* 15(b) tells us that the external force ("natural”) 
is not so far away from a sinusoidal form* 
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